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Summary 

This paper presents an analysis and sample re- 
sults of the lateral buckling and vibration of a com- 
pressively loaded column whose cross section is piece- 
wise constant along its length. The column is 
symmetric about its midspan and consists of three 
sections, with the center section having a stiffer cross 
section than the two identical outboard sections. 
Buckling and vibration characteristics of the column 
are determined from a numerical solution of the exact 
eigenvalue problems. Parametric structural efficiency 
analyses are performed using a nondimensionalized 
set of governing equations to determine the optimum 
ratio between the lengths of the center section and 
the outboard sections based on both buckling load 
and vibration frequency requirements. In these anal- 
yses two relationships between cross-sectional mass 
and bending stiffness are considered; one is a low- 
efficiency method for increasing the bending stiffness 
of the cross section, and the other is a high-efficiency 
method. The effect of axial load on vibration fre- 
quency is also examined and compared with that of 
a uniform column. 


Introduction 

Simply supported columns undergoing lateral vi- 
bration and Euler buckling have similar bending- 
moment distributions with low values near the ends 
of the column. This suggests that a column with a 
tapered cross section that is stiffer in the middle than 
at the ends can be more efficient (i.e., have less mass 
for a given buckling load) than a uniform column. A 
buckling analysis of a midspan, symmetric, uniformly 
tapered column is presented in reference 1 , and buck- 
ling test results for uniformly tapered columns appli- 
cable to large erectable space structures are presented 
in reference 2. These results have verified the in- 
creased structural efficiency of tapered columns com- 
pared with uniform columns. 

A disadvantage of the continuously tapered col- 
umn is that it is more difficult to fabricate than a uni- 
form column. The midspan symmetric column with 
a piecewise constant cross section shown in figure 1 is 
an approximation of a tapered column. This geom- 
etry potentially provides some of the increased effi- 
ciency of a tapered column without greatly increased 
fabrication complexity. This paper summarizes the 
results of a study to quantify the structural efficiency 
of a midspan symmetric column with a piecewise 
constant cross section based on both buckling load 
and vibration frequency considerations. The results 
derived herein are based on the solution of a non- 
dimensionalized set of governing equations and, 


therefore, are generally applicable to any column of 
this type. . 
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Figure 1. Simply supported column with piecewise constant 
cross section. 


Symbols 

Aj cross-sectional area of jth section 

Cjk constants in displacement solution for 

jth section 

(EI)j flexural stiffness of jth section 
1 1 length of outboard sect ions 

Z 2 half-length of center section 

M normalized mass 

P compressive load 

~P normalized compressive load 

P cr buckling load 

P CT normalized buckling load 

r radius of thin-walled circular cross 

section 


t 

Wj 

a 

0 

7 

Ajfc 


Pj 

r 

^ jk 
uj 


time 

lateral displacement of jth section 

longitudinal position coordinates 

ratio of length of center section to 
total column length 

ratio of flexural stiffness of center 
section to that of outboard sections 

ratio of mass per unit length of center 
section to that of outboard sections 

magnitude of fcth root of characteristic 
polynomial for jth section 

mass density of jth section 

thickness of thin- walled circular cross 
section 

fcth root of characteristic polynomial 
for jth section 

vibration frequency 



u) normalized vibration frequency 

Q 0 normalized vibration frequency of 

column with no axial load 

Derivation of Governing Equations 

For the present study it is assumed that the center 
section of the column shown in figure 1 is of length 
2/ 2 with flexural stiffness (£7)2, area A 2 , and density 
P2\ the outer sections are each of length l\ with 
flexural stiffness (El) 1 , area A\, and density p\. 
Because of symmetry in geometry and loading, it 
is possible to consider only one-half the column for 
analysis. Further, since the fundamental vibration 
and buckling modes of the column are expected to 
be symmetric, symmetry boundary conditions are 
assumed at the column midspan. Figure 2 shows 
this analysis model along with the assumed boundary 
conditions. 



An exact solution of this problem requires the 
solution of two fourth-order differential equations, 
one governing section 1 (0 < x < l\) and the other 
governing section 2 (l\ < x < l\ + l 2 ). In addition 
to the four boundary conditions given in figure 2, 
four more conditions arise from the enforcement of 
displacement, slope, moment, and shear continuity at 
the interface between sections 1 and 2 (x = l\). For 
convenience, a second coordinate axis (£) is located 
with its origin at the column midspan, and the 
displacement in section 2 is described as a function of 
this coordinate. The governing differential equation 
in the j th section is 

d*Wj P d 2 Wj PjAj d 2 Wj 
Ih* + (EI)j &r 2 + XElYj~W = ° 

The left-end boundary conditions in section 1 are 

w m = ^(o) = o ( 2 ) 


Similarly, the right-end boundary conditions in sec- 
tion 2 are 

( 3 ) 

Finally, the continuity conditions between sections 1 
and 2 are 


wi(Ji) - W2(h) 

(4) 

dw in ^ dw2,, , 

dx (il) = {h) 

(5) 

Mi 9 ^ 1 (k) = (Elh^h) 

(6) 


(7) 


Notice that the minus signs exist in equations (5) 
and (7) because the x and £ coordinate axes are in 
opposite directions. 

Equation (1) is a homogeneous partial differen- 
tial equation with constant coefficients. Assum- 
ing harmonic motion and applying separation of 
variables gives the following form for the lateral 
displacement: 


Wj(x,t) = Ce i “ t e' 1> i x (8) 

where C is an arbitrary constant. Substituting equa- 
tion (8) into equation (1) gives the following charac- 
teristic polynomial in tpj : 



(EI)j 


= 0 


( 9 ) 


so that 


fyl = ±i 
^2 = ± 


N 


P+ {P 2 + 4 Pj A j (EI) j cj 2 ) l/2 
2 (EI)j 




-P+ (P 2 + 4p ? Vl 7 (g/) 7 a;2) 1/2 


2 (EI)j 


( 10 ) 


(11) 


The right-hand side of equation (10) yields two 
purely imaginary conjugate roots. The right-hand 
side of equation (11) yields either a pair of real roots 
having the same magnitude and opposite signs or, in 
the static case (u> = 0), a pair of repeated roots equal 
to zero. 


Solutions 

The equations of the previous section can be 
solved to determine the critical buckling load of 
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the column by assuming the vibration frequency to 
be zero, or they can be solved to determine the 
fundamental vibration frequency for a fixed value of 
axial load. First, the solution of the buckling problem 
will be presented along with numerical examples to 
assess the structural efficiency of this type of column. 
Then, the solution of the vibration problem will be 
presented along with similar numerical examples for 
the case of no applied axial load as well as a range of 
axial loads. 

Solution for Column Buckling 

For u = 0, tp\2 = ^22 = 0) and the general 
solution can be written as follows for the two sections 
of the column: 


wi(x) = C\\ sin(Ajix) + C \2 cos(Anx) + C\$x + C\$ (0 < x < l\) 

( 12 ) 

W 2 (0 = C21 sin(A 2 if) + C22 cos(A 2 i£) + + C24 (0 < ^ < *2) 

(13) 

where Cn, C\ 2, C21, and C22 are arbitrary constants 
and An and ^21 are the magnitudes of the roots 
and ^21 of the characteristic polynomial given, 
respectively, by 


and 



(14a) 


(14b) 


Application of the boundary conditions in equa- 
tions (2) and (3) reduces the solutions to 


The existence of a nontrivial solution of equa- 
tions (17)-(20) requires that 


sin(An/i) sin(A2if 2 ) 


- cos(Anli) cos(A2i^2) = 0 
A 2 i 

( 21 ) 


which agrees with the result presented in refer- 
ence 1. Substituting the definitions for An and A 2 i 
from equations (14) into equation (21) gives the fol- 
lowing transcendental equation for the buckling load 


sin 


'Peril 

(El) 1 


1/2 

sin 


Pi 3-^2 

(El) 2 . 


(El) 2 


(El) i 


1/2 


[Ml 1/2 ( 

[(EI)i\ 


Peril 


(El) 2 


1/2 


= 0 ( 22 ) 


We now define the dimensionless parameters a, 0, 
and Per, respectively, as follows: 


a = 


13 = 


h 

h + h 

(El ) 2 

(EI)l 


Per = 


TT 2 (EI)l/4(l\ + h) 2 


(23) 

(24) 

(25) 


Notice that a is the ratio of the length of the center 
section to the total column length, 0 is the ratio of 
the bending stiffness of the center section to that of 
the outboard sections, and Pa- is the buckling load 
of the complete column normalized to that of a col- 
umn having the same length and a uniform bending 
stiffness equal to that of the outboard sections. Us- 
ing these nondimensional parameters, equation (22) 
can be rewritten in the form 


toi(x) = Cn sin (Ana:) + C131 (0 < x < l\) ( 15 ) 

w 2(£) = C22 cos (A21O + C24 ( 0 <£< (2) ( 16 ) 

Substituting equations (15) and (16) into the 
four continuity conditions (eqs. (4)-(7)) gives the 
following system of equations for C\\ , C13, C 22 , and 
C24: 

C\\ sin(Anii) + C13Z1 = C22Cos(A2ii2) + C24 (17) 

CnAn cos(Anii) 4- Ui3 = C22A21 sin(A 2 i/2) (18) 

(E/)iCnAi! sin(Auli) = (EPiCii^h cos^ik) (19) 
(fJ/JiCiiAjj cos(Anii) = (EI)i.C22^\ sin(A2i22) (20) 



Numerical Examples and Illustration of 
Mass Savings 

Because of its transcendental form, equation (26), 
with a and 0 prescribed, must be solved for Per nu- 
merically. The secant method (ref. 3) is selected be- 
cause of its simplicity and stable convergence char- 
acteristics. To ensure that the lowest root is found, 
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a search is performed to determine an initial interval 
in which the solution exists. During this search, the 
left-hand side of equation (26) is evaluated for succes- 
sively larger values of Per (starting with the lowest 
possible value, 1.0) until a sign change is detected. 
The initial solution interval is then bounded by the 
last two values selected for Further iteration 

yields the solution to the desired accuracy. 

In order to evaluate the structural efficiency of the 
column, it is necessary to calculate not only its nor- 
malized buckling load but also its normalized mass. 
For this, two additional dimensionless parameters 7 
and M are defined, respectively, as 


- P2 ^ 2 

piM 


(27) 


which is the ratio of mass per unit length of the 
center section of the column to that of the outboard 
sections, and 


— _ piAjh + P2Mh 

P\M{h + 12 ) 


1 + a ( 7 - 1) 


(28) 


which is the mass of the column normalized to that of 
a column having the same total length and a uniform 
cross section equal to that of the outboard sections. 

To illustrate the mass savings possible with this 
type of column, two different relationships between 
mass and bending stiffness are considered. The 
first is a low-efficiency method for increasing the 
bending stiffness of the cross section; the second 
is a high-efficiency method. In both cases it is 
assumed for simplicity that a single material is used 
throughout the column. In the first case the mass 
and bending stiffness are assumed to be proportional; 
thus 7 = (3 in equation (28). To illustrate this 
case, equations (29) and (30) give the approximate 
relationships for the area and moment of inertia of a 
thin-walled cylinder, where r is the radius and r is 
the thickness: 


Area « 2titt (29) 

Moment of inertia ~ 7 rr 3 r (30) 

Increasing the thickness while holding the radius 
constant increases the area and moment of inertia 
proportionally; thus 7 = (3. 

Figure 3 presents plots of the normalized col- 
umn mass versus normalized buckling load deter- 
mined from equation (26) for an array of length ra- 
tios (0.2 < a < 0.6 for fig. 3(a) and 0.6 < a < 1.0 for 
fig. 3(b)) and bending stiffness ratios (1.0 < (3 < 6.0). 
The dotted-line curves in these figures show the in- 
crease in mass necessitated by the increase in buck- 
ling load of a uniform column ( a = 1.0). Therefore, 


points that lie below these dotted-line curves repre- 
sent columns that are more efficient (lower in mass) 
than the uniform column. The data are separated 
into two parts for clarity. Notice that every curve of 
constant a is vertically asymptotic to the buckling 
load of a uniform column with a length 1 — a times 
the length of the original column, i.e., 2l\. This is 
consistent with the fact that when the bending stiff- 
ness of the center section is very large, it behaves as 
a rigid section. 


Uniform 0.5 



(a) 0.2 < a < 0.6. 



(b) 0.6 < a < 1.0. 

Figure 3. Mass versus buckling load where mass is propor- 
tional to bending stiffness. 
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For increases in buckling load up to a factor of 3, 
it appears that the optimum value for the length ra- 
tio a is approximately 0.7 because this corresponds 
to the lowest curve on the plot. Mass savings is 
calculated by determining the percent difference be- 
tween the lowest curve and the uniform column curve 
(ct = 1.0) at any given buckling load. It is seen 
that a mass savings of 10 to 12 percent can be real- 
ized by increasing the bending stiffness of the center 
70 percent of a column rather than the entire length 
of the column. Again, this conclusion is based on 
the assumption that mass and bending stiffness are 
proportional. 

In the second relationship considered between 
mass and bending stiffness, the mass is assumed 


to be proportional to the cube root of the bending 
stiffness, i.e., 7 = V/3 in equation (28). This 
case is representative of high-efficiency methods for 
increasing the bending stiffness of the cross section, 
and it is illustrated by a thin-walled circular cross 
section in which the radius is increased while the 
thickness is held constant. (See eqs. (29) and (30).) 

Figure 4 presents plots of the normalized column 
mass versus normalized buckling load for this case. 
As before, the dotted-line curves in this figure show 
the increase in mass using a uniform column (a = 
1.0), and the data points below these curves represent 
columns that are more efficient (less massive) than 
the uniform column. 



l 


0 1.0 2.0 3.0 4.0 5.0 

Normalized buckling load 

(a) 0.3 < a < 0.7. 

Figure 4. Mass versus buckling load where mass 



Normalized buckling load 
(b) 0.7 < a < 1.0. 

is proportional to cube root of bending stiffness. 


As was determined in the first case, it appears that the optimum value for the length ratio a is approximately 
0 7 In this case a mass savings of about 6 to 8 percent is realized by increasing the bending stiffness of the 
center 70 percent of a column rather than the entire length of the column. Again, this conclusion is based on 
the assumption that mass is proportional to the cube root of the bending stiffness. 


Solution for Column Vibration 


For vibration of the column under a specified axial load P, the general solutions are 

as follows: 


wi{x,t) = e iu3t jCn sin(Anx) + C u cos{\ n x) + C 13 sinh(A 12 :r) + C 14 cosh(Ai 2 x)] 

(0 < x < h) 

(31) 

w 2 (£, t) = e^ 1 \c 2 1 sin(A 2 iO + C 22 cos(A 2 iO + C 23 sinh(A 22 0 + C 24 cosh(A 22 £)] 

(0 < 4 < h) 

(32) 
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where, as before, C\\, C\ 2 - C 2 \, C 22 ■ C 13, C23, C14, and C24 are arbitrary constants and An, A12, A21, and 
A22 are defined as 


An = 


A21 = 


A12 = 


A22 = 


N 


P+ [P 2 + 4pii4 1 (£J) 1 «2] 1/2 

2{EI) l 


P + [P 2 + 4p 2 A 2 (EI) 2 u 2 ] 1/2 
2(EI) 2 


N 


-P + [P 2 + ApxAxiEI)^ 2 ] 112 

2(E/)i 


N 


-P + [P 2 + 4^yl2(g/)2^ 2 3 1/2 


(33a) 


(33b) 


(34a) 


(34b) 


2 (£7)2 

Upon application of the boundary conditions in equations (2) and (3), the general solutions reduce to 

w\(x,t) = e lult [Cn sin(Anx) + C13 sinh(Ai2x)] (0 < x < l\) (35) 

^ 2 (C i) = [C 22 cos(A 2 iO + C 24 cosh(A 2 2 ^)] (0 < £ < i 2 ) (36) 

Substituting equations (35) and (36) into the continuity conditions (eqs. (4)-(7)) gives the following system 
of equations that are written here in matrix form: 


sin(AnZi) 


sinh(Ai2/i) 


cos(AnZi) ^cosh(Ai 2 Zi) 


— COs(A 21?2) 
'Af[ sin (^21^2) 


-cosh(A22^2) 
^ sinh(A 2 2 ^) 


-sin(AnZi) (^) 2 sinh(A 12 /i) P cos(\ 2l l 2 ) -0 cosh(\ 22 l 2 ) 


-cos(An/i) 


cosh(Ai2/i) 0 


sin(A 2 i/2) P 


sinh(A22^2) 


'C n 
C 13 

C22 

(C24 


> = 0 


(37) 


where j3 is the bending stiffness ratio defined in 
equation (24). 

As in the case of the buckling problem, it is of in- 
terest to nondimensionalize the vibration eigenvalue 
problem to allow general solution curves to be con- 
structed. In addition to the dimensionless length, 
bending stiffness, and mass-per-unit-length parame- 
ters defined in equations (23)-(25), a dimensionless 
vibration frequency Q and compressive load P are 
defined, respectively, by 


OJ 


uj = 


7T 


jgfli 


P = 


1 tyMh + h) A 

p 

7 ^W)UW+W 


(38) 


(39) 


Hence, the vibration frequency and compressive load 
are normalized to the fundamental frequency and 
buckling load, respectively, of a uniform column of 
the same length with bending stiffness and mass per 
unit length equal to that of the outboard sections 
(ref. 4). 

All the parameters in equation (37) may now be 
written in_ terms of the five dimensionless quantities: 
a, /?, 7, P, and u). The resulting expressions are 


An/i = 


tt( 1 — a) 
2 



(40) 


^12*1 = 


7r(l — a) 
2 



(41) 
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^21*2 = 


A22^2 = 


not 

UP 


7 ra 



1/2 




1 / 2 ' 


1/2 


A 12 _ 

1 

1 

+ 

(?♦*: 

n 

An 

7 + ( 


1/2 


1/2 


(42) 


(43) 


(44) 


A21 _ J_ 

An y/fi 


P^(E1 + 7 /fo 2 ) 


7+lT 


1/2 


f + (^+“ 2 ) 


1/2 


1/2 


(45) 


A22 _ 

1 

1 

1 

4 - 

;? + ^) w ' 

An 


r 

+ 

1 


1/2 


(46) 


Numerical Examples for Column 
Vibration 

Equation (37) must be solved numerically to 
determine the normalized vibration frequency u>, 
given values for the other nondimensional param- 
eters. Again, the secant method was used to de- 
termine the lowest value of u) that causes the de- 
terminant of the matrix in equation (37) to vanish. 
The determinant of equation (37) was calculated us- 
ing a Gaussian elimination procedure presented in 
reference 5. 


the mass is proportional to the cube root of the bend- 
ing stiffness (i.e., 7 — 3 V^)- 

Figure 5 presents a plot of the normalized col- 
umn mass versus normalized vibration frequency for 
the first case (7 = /?). It should be noted that for 
a uniform column (a = 1.0), there is no change in 
the frequency as the bending stiffness and mass are 
increased proportionately. The significant result de- 
termined from figure 5 is that all column geometries 
considered exhibit frequencies that are lower than the 
uniform column having the same mass. Therefore, in 
the case where the distributed mass of the column is 
proportional to its bending stiffness, a column with 
piecewise constant cross section is less efficient than 
a uniform column, based on vibration frequency, in 
contrast to results based on buckling load. 



Figure 5. Mass versus vibration frequency where mass is propor- 
tional to bending stiffness. 


Natural vibration of columns without axial 
load . In the last section the structural efficiency of 
columns with piecewise constant cross sections was 
quantified based on buckling performance. It is also 
of interest to determine the vibration performance 
of these columns with no axial load. Consequently, 
equations (37) and (28) were solved to determine the 
vibration frequencies and masses for the same array 
of length ratios (0.2 < a < 1.0) and bending stiffness 
ratios (1.0 < 0 < 6.0) considered in the buckling 
solutions. As before, two cases were considered in- 
volving the mass per unit length of the cross section. 
In the first case it was assumed that the mass and 
bending stiffness of the cross section are proportional 
(i.e., 7 = (3). In the second case it was assumed that 


Figure 6 presents plots of the normalized column 
mass versus normalized vibration frequency for the 
second case (7 = 3 \/£). The dotted lines in these 
figures show the increase in mass for a uniform col- 
umn (a = 1.0); thus the data points below these 
lines represent columns that are more efficient than 
the uniform column based on fundamental vibration 
frequency. Again, the data are separated into two 
parts for clarity. 

The trends displayed in figure 6 are similar to 
those presented in figures 3 and 4 and discussed in the 
buckling analysis of the last section. In this case the 
total column mass can be minimized by stiffening the 
center 70 to 80 percent of the column. Furthermore, 
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the use of this design results in a mass savings of 
approximately 10 to 15 percent relative to a uniform 
column with the same vibration frequency. 

Vibration of columns under axial load . The 
effect of axial load on the fundamental vibration fre- 
quency of a uniform column is analyzed in refer- 
ence 6. The normalized formula for calculating this 
frequency can be written as 



where Q 0 is the normalized natural frequency of 
the column (no axial load), P is the normalized 
buckling load of the column, and P is the normalized 
axial compressive load in the column. 



(a) 0.2 < a < 0.7. 

Figure 6. Mass versus vibration frequency where 


Although it appears impossible to establish equa- 
tion (47) analytically for a column with a piece- 
wise constant cross section under axial load, nu- 
merical examples should shed some light on the 
question of its applicability. Table I presents the 
values of length, bending stiffness, and distributed 
mass ratios for five selected column configurations 
along with the numerical solutions for their normal- 
ized buckling loads and natural vibration frequen- 
cies. The sixth row in table_l lists the values se- 
lected for the axial load (P/'P cr ). (Note that the 
negative numbers indicate tension and the positive 
numbers indicate compression.) Finally, the last two 
rows of table I present the numerical solution for vi- 
bration frequency {(Q/Q 0 ) numerical) and the solution 
determined from equation (47) ((u>/u>< ? ) eq ( 47 j) for 
each configuration. It can be seen that these so- 
lutions agree very well, the differences likely being 
attributable to numerical inaccuracies. 



(b) 0.7 < a < 1.0. 

is proportional to cube root of bending stiffness. 


Figure 7 is a plot of equation (47) superimposed over the data points presented in table I. Although 
the set of configurations considered is certainly not exhaustive, it represents a reasonable range of the 
dimensionless parameters. The results strongly suggest that equation (47) indeed predicts the correct vibration 
frequency for columns with piecewise constant cross sections. If so, then to determine the effect of axial load 
on vibration frequency in a column of this type, it is necessary to calculate only the buckling load and natural 
vibration frequency of the column, then apply equation (47) to determine its vibration frequency under axial 
load. 
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Table I. Numerical Examples of Column Vibration With Axial Load 



Configuration 

Parameter 

i 

2 

3 

4 

5 

a 

0 

7 

Per 

&o 

P / P CT 

(^M>) numerical * * • 

(^M>)ea. (471 ' ' * 

0.4 

4.0 

4.0 

2.0107 

0.8463 

-0.75 

1.3167 

1.3229 

0.4 

4.0 

1.5874 

2.0107 

1.2354 

-0.5 

1.2225 

;.2247 

0.6 

2.0 

1.2599 

1.8123 

1.2162 

0.25 

0.8661 

0.8660 

0.8 

2.0 

2.0 

1.9738 

0.9970 

0.5 

0.7050 

0.7071 

0.8 

6.0 

1.8171 

5.5849 

1.7673 

0.75 

0.5004 

0.5000 



Figure 7. Effect of axial load on fundamental vibration 
frequency. 


Concluding Remarks 

The results of an analytical study of the buck- 
ling and vibration characteristics of a column with 
piecewise constant cross section have been presented. 
Parametric structural efficiency analyses determined 
that for increased buckling resistance, the optimum 
ratio between lengths of the stiffened center section 
and the entire column is approximately 0.7. Also, 
it was determined that a column using this ratio of 
lengths offers a mass savings of 6 to 12 percent rel- 
ative to a uniform column having the same buckling 


load. Furthermore, the magnitude of this mass sav- 
ings was shown to be dependent on the relationship 
between bending stiffness and mass per unit length 
of the column cross section. 

Similar parametric structural efficiency analyses 
were performed using a nondimensionalized set of 
the governing vibration equations. From these anal- 
yses it was determined that the relationship between 
bending stiffness and mass per unit length of the col- 
umn cross section has a great effect on the efficiency 
of the column from a vibration standpoint. If the 
mass per unit length and bending stiffness are pro- 
portional, a column with a piecewise constant cross 
section is less efficient than a uniform column, based 
on fundamental vibration frequency. However, if the 
mass per unit length is proportional to the cube root 
of the bending stiffness, a column with a piecewise 
constant cross section is more efficient than a uni- 
form column, based on either fundamental vibration 
frequency or buckling load. 

Finally, numerical results strongly suggest that 
the relation between axial load and fundamental 
vibration frequency for a uniform column also holds 
for a column with piecewise constant cross section. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
January 31, 1991 
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